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ABSTRACT 

We are conducting a search for supermassive black holes (SMBHs) with masses below 
~ 10^ Mq by looking for signs of extremely low-level nuclear activity in nearby galaxies 
that are not known to be AGNs. Our survey has the following characteristics: (a) X-ray 
selection using the Chandra X-ray Observatory, since x-rays are a ubiquitous feature 
of AGNs; (b) Emphasis on late-type spiral and dwarf galaxies, as the galaxies most 
likely to have low- mass SMBHs; (c) Use of multiwavelength data to verify the source 
is an AGN; and (d) Use of the highest angular resolution available for observations in 
x-rays and other bands, to separate nuclear from off- nuclear sources and to minimize 
contamination by host galaxy light. Here we show the feasibility of this technique to 
find AGNs by applying it to six nearby, face-on spiral galaxies (NGC 3169, NGC 3184, 
NGC 4102, NGC 4647, NGC 4713, NGC 5457) for which data already exist in the 
Chandra archive. All six show nuclear x-ray sources. The data as they exist at present 
are ambiguous regarding the nature of the nuclear x-ray sources in NGC 4713 and NGC 
4647. We conclude, in accord with previous studies, that NGC 3169 and NGC 4102 are 
almost certainly AGNs. Most interestingly, a strong argument can be made that NGC 
3184 and NGC 5457, both of type Scd, host AGNs. 
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Introduction 



The past decade has seen extraordinary improvement in our understanding of supermassive 
black holes (SMBHs) — their growth and evolution, and their links with their host galaxies. We 
now realize that galaxies hosting SMBHs at their centers are the rule rather than the exception. 
The observed correlations of the masses M, of the SM BHs with properties of th e ir host galaxies 



for e xample with the bulge stellar velocity dispersion (iFerrarese &: Merritt 2000 : Gebhardt et al 



2000), show that there is a close link between the formation and evolution of galaxies and of the 
SMBHs they host. Furthermore, a comparison of the SMBH mass function required to explain the 
observed luminosit y function of active galactic nuclei (AGN s) with estimates of the local SMBH 
mass function (e.g. iMarconi et al.ll2004l : IShankar et al.ll2004l ) shows that not only must SMBHs be 
very common in massive galaxies but that most, if not all, of these black holes are relics of AGNs 
active in previous epochs. Thus knowledge of the local SMBH mass function enables us to put 
constraints on theories of galaxy and SMBH formation and growth, and AGN lifetimes. 

Estimates of the local SMBH mass function are anchored by reso lved stellar or gas dyna mical 
measurements of the masses of ~ 30 SMBHs (see, e.g., the review by IFerrarese fc: Fordll2005l ). and 
otherwise based on the distribution of host galaxy properties (luminosity of the bulge or bulge 
stellar velocity dispersion a) and known scaling relationships between t he mass of the SMBH and 
these properties (most prominently M, — a and M, — Mbuige ; but see I Greene &: Hoi l2007al for a 
more direct estimate ). Most of the measured SMBH masses are ~ 10^ Mq or greater, however, 
as the sphere of influence of a less massive SMBH is extremely hard to resolve even at moderate 
distances, even with HST. For example, the sphere of influence of a 10^ Mq SMBH at 15 Mpc is 
~ 30 milliarcseconds (mas). As a result, while different estimates of the mass function agree for 
lO^ MfTi < M, < lO^MfTj , the low-mass end (Af, < 10^ M0) often has discrepancies (see, e.g.. Fig. 7 
of iGraham et al.ll2007l for a comparison of different authors' estimates of the mass function). A 
second source of uncertainty at the low-mass end is the fact that it is unknown how the scaling 
relationships extrapolate to very late-type spiral galaxies, which have little or no bulge component, 
and to very low mass galaxies (dE and dSph). Yet SMBHs do exist in very late-type sp irals, e.g. 
NGC 4395, a spiral galaxy of type Sdm, with M, ~ 3 x lO^M© dPeterson et al.l boosl'l. and in 
very low mass galaxies, e.g. POX 52, a dwarf galaxy, with M, ~ 3 x 10^ Mq (jBarth et al.l 12004 ). 
Questions that naturally arise at this point are: Do the scaling relationships break down at low 
masses? What determines the mass of an SMBH: the mass of the bulge or the mass of the dark 
matter halo? Is there a lower bound to the local SMBH mass function? A well defined sample of 
low-mass SMBHs is needed to answer these questions. 

Since we cannot detect low-mass SMBHs by their dynamical signature, looking for them by 
signs of their accretion activity may be the only viable way of detecting them. This of course 
limits detection to the subset of low-mass black holes that are active, but this fraction can be 
expected to be large, for the following reason. We now understand that the "quasar era" is a func- 
tion of luminosity, with the space density of the most luminous quasars peakin g at high redshift 
and that of lower luminosity quasars peaking at progressively lower redshifts (jFiore et al.ll2003l : 
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Hasinger et al.ll2005l ). This is often referred to as the "downs izing" of AGN activ i ty with cosmic 
epoch. Moreover, at least some models of black hole growth ([Marconi et al.ll2004l : lMerlonill2004l ) 
require anti-hierarchical growth. That is, higher mass SMBHs attain most of their mass at high 
redshift while lower mass SMBHs grow at progressively lower redshifts. The trends of AGN down- 
sizing and anti-hierarchical growth, extended logically to the smallest mass SMBHs, imply that 
these objects were active in recent times and may still be accreting at the present epoch. 

The Eddington luminosity of a IO^Mq SMBH is only ~ 10^^ erg s~^; low-mass SMBHs, even 
if accreting at high rates, will be low-luminosity AGNs (LLAGNs). The converse is not true; that 
is, not all LLAGNs have a low-mass SMBH. In addition to a low SMBH mass, the low luminosity 
of an LLAG N may be caused by a low rate or radiatively inefficient mode of accretion to an SMBH 



of any mass (jSoria et al 



2nn6. 



present a study of massive SMBHs accreting at very low rates). 
Finally, obscuration may further lower the observed luminosity. Low-ioniza tion nuclear emission- 
line region (LINER) nuclei have been studied at multiple wavelengths (e.g. Eracleous et al. 200^; 



Satyapal et al. 



2004 



Dudik et al. 



2005 



Flohic et al.l 120061 ) to identify LLAGNs among them, and 



these studies have demonstrated that these AGNs are not necessarily the same type of object with 
the same physical characteristics. Since we are interested specifically in the low-mass end of the 
SMBH mass function, it is necessary to identify among the AGNs those that can be expected to 
have the smallest black holes. One approach is to use mass estimators based on the luminosity 
of the AGN and the width of the broad compon ent of emission lines in its spectrum. This is the 
technique that was used bv lGreene &: Hoi (|2007bl ) in constructing their sample of low-mass SMBHs. 
A second approach, and the one we use, is to look for AGN activity in galaxies (late-type spirals, 
dwarf galaxies) where the known host galaxy-SMBH scaling relationships predic t the lowest-mass 
SMBHs would reside. This approach has been used by lSatyapal et al.l (|2007l . 120081 ) to find candidate 
low-mass SMBHs. However, the scaling relationships are only statistical and cannot be used to 
estimate the mass of a particular SMBH; thus the second approach requires an independent estimate 
of the SMBH mass. 

A search for active low-mass SMBHs requires confirmation of the presence of an AGN in each 
candidate nucleus, and measurement of the mass of the SMBHs in the confirmed AGNs. The 



method used bv lGreene &: Hoi (l2007bl ) has the virtue of effectively combining all of the above into 
a single step. However, as the luminosity of an AGN decreases, the optical spectrum of the galaxy 
nucleus becomes more and more dominated by host galaxy light, and the signature of the AGN 
becomes difficult to detect. Even when the optical spectrum shows no clear evidence of an AGN, 
however, such ev idence may still b e present in other wavelengths, such as x-r ay an d radio (e.g. 
Filho et al.l[2004l l. and infrared (e.g. bale et"al]l20od : Isatvapal et ahlbooj boosl l (See Eol[2008l for 
a review of nuclear activity in nearby galaxies). For the lowest-luminosity AGNs, therefore, it is 
possible that a system based on optical spectra would not classify the nuclei as AGNs at all. We 
choose to use x-ray selection to identify candidate AGNs for the following reasons: First, x-rays 
can penetrate obscuring material which may be hiding the line emitting regions. Second, there are 
fewer sources of x-rays in a galaxy than there are of optical and UV emission and so dilution of the 
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AGN signature by host galaxy light is less of a problem. Where the luminosity of the AGN is low 
to begin with, even a modest amount of absorption may result in the signal being b elow the optical 



back g round iniposed by the galaxy. Third, even if, as expected in some theories (jEracleous et al 



I995I : iNicastrd l200d : iNicastro et al.l 120031 : iLaorl 120031 ) AGNs that have luminosities or accretion 
rates below a cut-off value do not have broad-line regions, they should still be detectable in x-rays. 
X-ray observations have in fact detected AGNs in what were thought to be "normal" galaxies (e.g. 
Martini et al.ll2002l ). The disadvantage, as discussed in ^ is that x-ray observations by themselves 
cannot always distinguish between AGNs and other x-ray sources, such as x-ray binaries (XRBs) 
and ultraluminous x-ray sources (ULXs). Multi- wavelength data are needed to determine the type 
of source. 

As the first step towards assembling a sample of low-mass SMBHs, we are conducting an x-ray 
survey of nearby (within 20 Mpc), quiescent spiral galaxies using the Chandra X-ray Observatory 
to look for nuclear x-ray sources, with an emphasis on late-type spirals. The high angular resolution 
of Chandra is necessary to ensure that any detected source is really at the center and is not an off- 
nuclear source. The survey is sensitive to an (unobscured) SMBH of mass M, = W^Mq radiating 
at ~ 2 X 10^^ L^dd out to the survey limit. We will present details of the survey sample selection and 
the Chandra observations in a future paper. Here, we present the methods used and the feasibility 
of detecting AGNs using these methods by applying them to six galaxies that meet selection criteria 
similar to those used for the survey sample, and for which x-ray data already exist in the Chandra 
archive. 

The paper is organized as follows: ^ describes the criteria used to select the six galaxies 
presented here; ^ describes the observations and data analysis, with individual targets discussed 
in S^SHHSISl the results are discussed in 331 



2. Sample selection 



Starting with Third Reference Catalogue of Bright Galaxie^ (RC3; de Vaucouleurs et al 



I991I ). we imposed the following main filters. (1) Morphological type: we selected spiral galax- 
ies of type Sa through Sd (1.0 < T < 7.0); (2) Distance: galaxies that had recession velocity 
cz < 3000 km s~^; (3) Galactic latitude: |6| > 30° to avoid x-ray absorption by gas in our own 
Galaxy; (4) Inclination: we required log{a/b) < 0.4, where a and b are the projected lengths of 
the semi-major and semi-minor axes, respectively. This was to avoid obscuration by the disk of 
the host galaxy; (5) Nuclear inactivity: the galaxy must not be known to be an AGN or star- 
burst. Starbursts were excluded since the vigorous star formation also increases the likelihood of 
the presence of XRBs and ULXs. LINERs were included as there is still debate about whether the 
fundamental source of energy in these nuclei is an AGN or star formation. We also required that 



In practice, we accessed the RC3 table in the SDSS database. 
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the object have the value of the "goodposition" flag set to 1 in the RC3 table in the Sloan Digital 
Sky Survey (SDSS) database, which indicated that the coordinates of the galaxy were accurate to 
0.1 s in RA and 1" in Dec. All objects passing these filters were cross-correlated with SDSS Data 
Release 3, with a match radius of 30". This resulted in a list of 76 galaxies. This list was then 
checked against publicly available data in the Chandra archive. The final sample consists of six 
galaxies: NGC 3169 (Sa), NGC 3184 (Scd), NGC 4102 (Sb), NGC 4647 (Sc), NGC 4713 (Sd), and 
NGC 5457 (Scd). 



3. Data analysis 

Details of the observations that were analyzed are given in Table [TJ X-ray data were down- 
loaded from the Chandra archive and analyzed using version 3.4 of Chandra Interactive Analysis 
of Observations (CIAO) software. In each case the Level 2 event list was processed using the 
observation-specific bad-pixel file and the latest calibration files, and filtered to exclude times when 
the instrument experienced a background fiare. The event lists were then converted to FITS images 
and the CIAO wavelet source detection tool wavdetect was run on them to determine source posi- 
tions. For each observation, source counts were extracted from a circle centered on the wavdetect 
source position, and with a radius equal to the greater of 4.67 pixels (2.3") and the 95% encircled- 
energy radius at 1.5 keV on ACIS-S at the position of the source. Background counts were taken 
from an annulus with an inner radius of twice, and an outer radius of five times, the source circle 
radius after excising any sources that happened to fall in the background region. Counts were 
extracted in the 0.3 - 8.0 keV (Broad), 2.5 - 8.0 keV (Hard) and 0.3 - 2.5 keV (Soft) bands. A 
hardness ratio (HR) was defined as HR = {H — S)/{H + 5), where H and S are the counts in 
the hard and soft bands, respectively. Spectral fits were performed with the CIAO tool Sherpa. 
Uncertainties reported in fit values represent the 90% confidence level for one parameter. In all 
estimations of the AGN Eddington ratio we have assumed that the bolometric luminosity is ten 
times the 2-10 keV luminosity of the AGN. 

Where available, optical and UV data from HST (WFPC2) were downloaded from the archive 
and analyzed using IRAF to extract nuclear fluxes. We note here that we use "nuclear flux" to 
refer to the total background-subtracted flux at the location of the nucleus. We use this flux, 
and colors if observations in more than one band exist, only as a consistency check, not to obtain 
photometry of any putative AGN. At these extremely low luminosities (Lboi ~ 10^^-10^*^ erg s""*^), 
the contribution of the AGN, even if present, will be a small fraction of the nuclear flux in the 
optical and UV bands. In most cases, therefore, we expect the fluxes and colors to be consistent 
with the absence of an AGN. We check for anomalous colors or fluxes, since these may indicate the 
presence of an additional source of flux, possibly an AGN. As such, we perform simple aperture 
photometry for resolved nuclear sources, with the aperture size matched to the source size. For point 
sources background-subtracted source counts were determined using imexam, using a 2-pixel radius 
aperture for the source, and an annulus with an inner radius of 7 pixels and an outer radius of 12 
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pixels for the background. These were corrected for CTE effe cts (lDolphinll2002ll the n extrapolated 
to a 0.5" radius aperture using the encircled-energy tables of lHoltzman et al.l (|l995l ). A correction 
of —0.10 mag was added to represent the magnitude in an infinite aperture. Proprietary ACS data 
for NGC 4713 were made available by Martini et al. (2008, in preparation). The flux of the resolved 
nuclear source was estimated using aperture photome try on the drizzled im age. We use infrared 
fluxes from the Two-Micron All Sky Survey (2MASS: ISkrutskie et al.l 120061 ) Point Source Catalog 
(PSC) for similar consistency checks. 

Individual targets are discussed below. For each target, we first describe the observed prop- 
erties in the x-ray, and then in other wavebands when such observations exist. We then consider 
whether the evidence supports the hypothesis that the x-ray source is an AGN. 



3.1. NGC 3169 



This is a galaxy of type Sa, at a distance of 19.7 Mpc (lTullylll988l ). and its nucleus is classified 
as a LINER (IHo et al.lll997al^. The Chandra observation of the nucleus of NGC 3169 was analyzed 
by iTerashima &: WilsonI (120031 ) and the nucleus identified as a low luminosity AGN, but the x-ray 
data are analyzed and presented here again for completeness. NGC 3169 was observed by Chandra 
for 2 ks and detected with 159 counts. Its hardness ratio HR = +0.86 makes it the hardest of the 
six sources discussed here and the only one to have a positive HR. Its spectrum, shown in Fig. [H 
can be reasonably fit by an absorbed power-law, with Nfj ~ 10^'^ cm~^ and F ~ 2. Spectral fit 
parameters are shown in Table [H with the first row showing a fit using the Cash statistic and data 
binned to five counts per bin, and the second row a fit using the Gehrels statistic and data 
binned by 20 PHA channels. Best-fit parameter values from both fits are con sistent within errors 
with each other and with the results obtained by ITerashima &: WilsonI (120031 ) . Fluxes used in the 



analysis are from the Cash fit. The unabsorbed broad band fiux is /(0.3 — 8keV) = 1 x 10 erg 



cm ^ s ^, implying a luminosity -L(0.3 



8keV) = 5 X 10^^ erg s ^ for the assumed distance. 



Observations with the Very Large B aseline Array (VLBA) show the nucleus is a mas-scale 



radio source at 5 GHz (iNagar et al.ll2005l ). The nuclear fiux density is /(5 GHz) = 6.6 x 10 erg 



cm s Hz , or uLi/ — 1.5 X 10'^'^ erg s ^. The m easured brigh tness temperature exceeds 10 



7.7 



K and rules out starbursts and supernova remnants (INagar et al.ll2005l ). 



The nucleus was detected by 2MASS. The observed magnitudes J = 11.2, H = 10.5, Ks = 10.0 
imply luminosities uL^, ~ (4 — 6) x 10^^ erg s~^ in these bands. 



3.1.1. The nature of the nuclear emission 



The combination of being a sub-parsec scale radio source as well as a luminous hard x-ray 
source points to the source being an AGN. Further evidence comes from the high luminosity. The 
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source is obscured, as made clear from the x-ray spectrum. The infrared luminosity therefore is a 
better indication of the true bolometric luminosity. This luminosity, a few xlO^^ erg s~^, is higher 
than is expected from XRBs or nuclear star formation regions. 



Dong &: De RobertisI (120061 ) perform a bulge-disk decomposition of NGC 3169 and use the 
2MASS K-band lumi n osity of the galaxy to estimate a central black hole mass log (M,/Mq) = 8.2. 
Heraudeau &: SimienI (jl998l ) report the central stellar velocity dispersion cr^, = 163 km s"^, which 
implies log (M./M©) « 7.8. Assuming log (M./Mo) ~ 8 for simplicity, the corresponding L^dd ~ 
10^^ erg s~^. The inferred Lboi ~ 10^^ erg s~^ then implies LboiZ-^^Edd ~ 10~^. 



3.2. NGC 3184 



This is of type Scd, classifi ed as having an H II nucleus by IHo et al.l (jl997al ) , and is at a 



distance of 8.7 Mpc (|Tullvlll988l ). It was observed with Chandra twice, one month apart, for 40 
ks and 25 ks. Both observations used the ACIS-S3 chip. The nuclear x-ray emission appears to 
consist of two components (Fig. [2]), a fainter, point-like source in the north and a brighter, more 
extended source in the south. The total extent of the nuclear emission is about 4" (~ 170 pc). 
Both components are soft and have HR ~ —0.8. 

Since the northern component has only ~ 30 counts, detailed spectral fitting is not possible. 
We check for consistency of the data with different models using unbinned data. An absorbed 
power-law can be fitted to the data, and despite the degeneracy in the parameters, an intrinsic 

^ 6 X 10^^ cm~^ is inconsistent at the 3-cr level. Fixing Nfj at the Galactic value gives 
a range of acceptable power-laws with 1.2 < F < 2.8 (3-cr limits). An unabsorbed blackbody 
with temperature 0.28 keV < kT < 0.56 keV (3-0" limits) is also consistent with the data. The 
observed broadband luminosity of this component, corrected for Galactic absorption, corresponds 
to -Z^o.s-SkeV ~ 2 X 10^'^ erg s~^. The difference in count rates between the two observations is not 
statistically significant, so this source does not vary on month timescales. 

There are archival HST observations in the UV (WFPC2, F300W) and optical (WFPC2, 
F606W) of this galaxy. The F606W image has the nucleus on the WF4 chip, and clearly shows a 
central point source. The nucleus is on the PC chip in the F300W image. There is no central point 
source detected. Both images also show diffuse nuclear emission and a prominent spiral arm in the 
central 1". Other than the nuclear source in the F606W image, there are no obvious counterparts 
of the x-ray sources in either the UV or the optical. We note here that the coordinates in the HST 
images are misaligned with each other and with the x-ray image. The images were therefore aligned 
by cross-correlating point sources in each individual image with SDSS sources that are classified as 
stars. The aligned HST images match each other to within 0.5", and match the Chandra image 
to within 0.7". These offsets are comparable to the absolute astrometric accuracy of both Chandra 
and HST. 
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3.2.1. The nature of the nuclear emission 



The fact that we see x-ray emission in the soft band but almost none in the hard band argues 
against the emission being the continuum from the AGN, whether direct or scattered. The source 
is not seen in the UV image. Both of these points suggest the presence of heavy obscuration. 
The observed x-rays could instead be soft emission from circumnuclear ionized gas, if the AGN is 
completely obscured and what is visi ble is mostly re-processed radiation. This is in fact the case i n 



some heavily obscured Seyfert 2s (e.g. iBianchi et al. 



2006 



Levenson et al.ll2006l : iGhosh et al.ll2007l ) 



where the dominant emission is the soft emission from the Narrow-Line Region. The Chandra and 
HST observations, therefore, are inconclusive regarding whether the source is an AGN, but do not 
rule out that possibility either. 

A stronger argument that the source is an AGN derives from infra-red data. T hi s galaxy is 



part o f the Spitzer Infrared Nearby Galaxies Survey (SINGS: iKennicutt et al.l 120031 ). iDale et al 



(|2006l ) have used the equivalent width of the PAH feature at 6.2/xm and the fluxes in a mix of 
high- and low-ionization lines ([S IV] 10.51/xm, [Ne II] 12.81/im, [Ne III] 15.56//m, [S III] 18.71//m, 
[O IV] 25.89//m, [S III] 33.48/im, [S II] 34.82^m) to create diagnostic diagrams that distinguish 
between AGN and star-forming galaxies. This nucleus falls into the "transition" region between 
AGNs and H II regions. This suggests an AGN component to the emission exists that may have 
been diluted because of the large aperture used (~ 20") to extract the fluxes. In IRAC images the 
nuclear source is resolved (~ 3.5" FWHM compared to the ~ 1.7" FWHM of the PSF). Nuclear 
fluxes were extracted using 3" apertures in the IRAC channels (D. A. Dale, private communication). 
Host galaxy emission appears to dominate the MIPS (24, 70, and 160 ^m) fluxes, but the observed 
IRAC colors, [3.6] - [4.5] = -0.26 ±0.16 and [5.8] - [8.0] = +0.59 ±0.16 (magnitude s in AB^ 
system), are redder than more than 80% of normal late- type galaxies (lAssef et al.ll2008l ). This is 
expected if there is an AGN, as AGN power-law emission falls off more slowly than galactic emission 
in the NIR. Thus, the IR line ratios and IRAC colors strongly argue for the source to be an AGN. 

The nucleus was also detected by 2MASS. The J, H, and Ks magnitudes, from the 2MASS 
Point Source Catalog (2MASS PSC), are J = 13.3, H = 12.7, and Ks = 12.5, which correspond 
to uLi^ ~ 10^^ erg s~^. Such a luminosity can be easily produced by an AGN. While all the 
observations are consistent with there being an AGN in NGC 3184, we cannot rule out a nuclear 
super-starcluster as the source. The x-ray emission could be from one or more XRBs in such a 
cluster. The x-ray luminosity (~ 10^^ erg s~^) is in the range seen from XRBs. The pres e nce of a 



nuclear star format i on reg ion is indicated by the optical line ratios observed by lHo et al.l ()1997al ). 
In addition, iLarsenl (|2004l ) reports a candidate nuclear star cluster in NGC 3184. 



In conclusion, there are two scenarios that are consistent with the observations. In both there 
is a nuclear star cluster that dominates the optical emission. In the first, there is no AGN and the x- 
rays are produced by one or more XRBs. In the second, there is a low-luminosity AGN in the center 
of the star cluster, so heavily obscured that we do not observe any direct or scattered emission. 
The x-ray emission arises from photoionized gas immediately surrounding the AGN. Using the 
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luminosity derived using the Bremsstrahlung model above, and assuming that approximately 1% 
of the AGN luminosity is reprocessed into the plasma x-ray emission, the AGN has a luminosity 
~ 10^^ erg s~^. The infrared line ratios and IRAC colors from Spitzer strongly argue for this 
scenario. 



3.3. NGC 4102 



This is an Sb galaxy, classified as H II in iHo et al.l (|l997al ) a nd as H II/ LINER in the 



NASA/IPAC Extragalactic Database (NE D), at a distance o f 17 Mpc (Itu11vIIi988|). The Chandra 



observ ation of this galaxy was presented bv lDudik et al.l (|2005l ) and also bv lTzanavaris Sz Georgantopoulos 



(|2007l ) . There are about 350 counts in the x-ray image of the nucleus. As a hardness ratio map 



(Fig. [3]) shows, the hard and soft emission are well segregated into a fairly hard (HR = — 0.20ib0.08), 
point-like source with extended, very soft (HR = -0.88l[J:[]^) emission to the west of it. Wavdetect, 
run on hard and soft band images, detects a small hard source and an encompassing soft source 
whose centers are 1" apart. We therefore fit the spectra of the harder "core" and the softer ex- 
tended emission separately, using counts from the regions shown in Fig. [3l The spectrum of the 
core has 171 counts, and shows a broad line between 6 and 7 keV, suggestive of reflection. The 
spectrum was fit after binning to 5 counts per bin, using the Cash statistic (cstat in Sherpa). A 
simple absorbed power-law model, fit to the energy range 0.3-3 keV, shows increasingly positive 
residuals at energies above 3 keV, again suggestive of a reflection component. The fit quality is 
poor, with a statistic value of 200 for 27 d-o-f. Using a reflection model (xspexrav) reduces the 
statistic to 60 for 26 d-o-f. The possible line at 6.4 keV is not well constrained because of the lack 
of counts on the high energy side of the line (there are 8 counts with E > 7 keV). For the line, we 
added a Gaussian at the fixed position of 6.4 keV and of fixed FWHM 0.3 keV, allowing only the 
amplitude to vary. This reduces the statistic further, to 40 for 25 d-o-f. The final parameters are 
r = 2.2+[j| and reflection scahng factor R = 129^8^. The formal equivalent width (EW) of the 
line is 2.5 keV, but the uncertainty in the amplitude of the line is of the order of ~ 50% and in 
the normalization of the reflected component, ~ 30%. The flux, excluding the line and corrected 
for Galactic absorption, is -F(0.3 — 8keV) = 4.2 x 10"^'^ erg cm~'^ s~^. A fit using the statistic 
(shown in Table H]) gives best-fit parameter values consistent with those from the Cash fit, with 
the exception of the equivalent width of the line, which is 3.3 keV in this fit. Given the poor 
determination of the continuum near the line the difference is not significant. The equivalent width 
of the Fe Ka in type 1 AGNs is typically of the order of 300 eV, and a large EW implies that the 
direct continuum is suppressed and the source is reflection dominated. In NGC 4102 the formal 
EW is an order of magnitude larger; therefore even after accounting for the large uncertainty we 
can conclude that the true EW is greater than the typical type 1 value. 

The extended emission is very soft, with just 5 counts (out of 115) above 2.5 keV. The spec- 
trum was fit after binning to 5 counts per bin and using both the and Cash statistics. A 
bremsstrahlung model provides an acceptable fit (x^ = 9.6 and Cash statistic = 30 for 18 d-o-f), 
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and the MEKAL model is less favored (x^ = 15.6 and Cash statistic = 41 for 18 d-o-f). Best-fit 
parameter values obtained using the two statistics are consistent within the errors. The best-fit 
temperature kT ~ 1 keV is higher than what may be expected of 100 pc-scale circumnuclear gas. 
The plasma surrounding an AGN may be photoionized rather than collisionally ionized, and there 
may be complexity in the spectrum that is hidden because of the poor quality. The flux in the ex- 
tended emission is F(0.3— 8keV) ~ 2 x 10^^'^ erg cm~^ s~^. Spectral models and best-fit parameter 
values for both the core and the extended emission are shown in Table [H The fits for both the 
core and extended components are shown in Fig. [H 

HST imaging observations in visual and near-infrared (NIR) bands show that the nucleus is 



clum py and has large amounts of dust, and that there is circumnuclear star formation (| Car olio et al 



19971 ). High absorption and reddenii ig towards the n ucleus is also indicated by the Balmer decre- 
ment which gives E{B — V) = 1.00 (|Ho et al.lll997al ). The nucleus was detected by 2MASS. The 
2MASS PSC gives J = 10.9, H = 9.8, and = 9.2. For a distance of 17 Mpc, these correspond to 
vL^{J) w 6x10^2 ergs-\ uLy{H) w 8x10"^^ ergs-\ &nduLy{Ks) ^ 7x10^^ erg s"^ The nucleus 
is highly luminous in the far- in frared (FIR). Infrared As tronomical Satellite (IRAS) observations 



indicate -LpiR ~ 10 erg s (|Mouri Taniguchil Il992l ). where LpiR is the luminosity between 



40/xm and 120/xm, albeit in the arcminute-scale IRAS apertures. The nucleus was also detected 
in radio by the FIRST survey (1.4 GHz). The stated radio position is within 0.5" of the Chandra 
position, and therefore within Chandra astrometric accuracy. The radio beam size, however, was 
3.75" X 2.84". The integrated radio flux density was 223 mJy, with peak flux densi ty of 167 mJy. 
This is the same integrated 1.4 GHz flux density r e porte d by ICondon et al.l (|l982l ). fifteen years 
prior to the 1997 FIRST observation. I Condon et al.l (119821 ) also report a higher resolution 4.9 GHz 
observation. The 4.9 GHz map shows radio emission extended in the northeast-southwest direction 
with a total extent of about 5.2" (~420 pc), and with two peaks separated by ~0.9" (~70 pc) and 
also aligned along a NE-SW axis. The position of the southwest peak coincides with the hard x-ray 
"core" seen in the Chandra observation. The radio spectral index a = —0.7 {Su cx i^°) suggests a 
synchrotron rather than thermal origin for the radio emission. 



3.3.1. The nature of the nuclear emission 



We first consider the case where there is no AGN and all of the observed emission is due to 
star formation, and in particular consider the radio emitting region, which, as mentioned above, 
has a size of the o rder of 400 pc. If this region obeys the radio-FIR correlation for non-AGN galax- 



ies ( Condonlll992l ). then it has an FIR luminosity of ~ 10^^ erg s ^, exactly what was measured 



by IRAS. Therefore the radio and FIR observations are consistent with each other and with the 
star-formation-only hypothesis, requiring only t he assumption that this region dominated the FIR 
emission within the IRAS aperture. Following ICondonI (119921 ) . we use the radio emission to esti- 
mate a supernova rate fsN and consequently a star formation rate (SFR). The luminosity density 
measured by FIRST, L;^(1.4GHz) 8 x 10^* erg cm~^ s~^ Hz~^, implies fsN ~ 0.08 yr~^, and 
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SFR(M > 5Mq) 2Mq yr-^ In the samples of lGrimm et alJ ((20031) and lRanalli et alJ ^2Q0i ). the 
galaxies that have this star formation rate have x-ray luminosities ranging from a few xlO^^ erg 
s~^ to a few xlO^'^ erg s~^. The observed nuclear x-ray luminosity of NGC 4102, Lx ~ 1-5 x lO^'^ 
erg s~^ is really a lower limit because the amount of absorption is undetermined, but for the pur- 
pose of the a rgun ient here may be considered to be consistent with the Lx~SFR relationship of 
Ranalli et alJ and iGrimm et alJ On the face of it, therefore, all of the observations are consistent 
with a starburst origin to the emission, but this explanation requires an extraordinarily intense 
starburst in a non-interacting galaxy. The radio, FIR, and x-ray luminosities require 2 Mq yr~^ 
of massive (M > 5Mq) star formation within a 400 pc region. For comparison, the massive star 
forma tion rate in the interac ting starburst galaxy M 82 integrated over the whole galaxy is 2.2 Mq 
yr-i dNeff Ulvestadlboool ). and the total SF R (M > 5M(t,) in the merging pair NGC 4 038/4039 



is estimated to be between 5 and 10 Mq yr ^ (jNeff &i Ulvestad 



2000 



Grimm et al.ll2003l ^. 



In favor of the AGN hypothesis, NGC 4102 shows an approximately conical region of outflowin g 
gas that has a higher [O IIl]/H/3 ratio than the surrounding star forming regions (jCanda et al.ll2006l ). 
suggesting exposure to a harder ionizing radiation than stel lar continuu m emission, and reminiscent 
of the ionization cones sometimes seen in Seyfert 2s (e.g. |Poggdll988l ). The [O III] line profile in 
this region is broa de r than in the su r round ing regions and cannot be fit with a single Gaussian 
(jCanda et al.l 12000 ). iGongalves et al.l (jl999l ) reported a broad component (FWHM f« 560 km s~^) 
to the [O III] line as well. The broad component is weak, comprising just 5-7% of the flux in the line 
in the case of the Ha and II/3 li nes. Using this co mponent gives [N II]A 6583/IIa = 1.57, similar to 
LINERs and Seyferts, and thus iGoncalves et al.l argue for the presen ce of a very w e ak Se yfert 2 in 
NGC 4102. In this respect NGC 4102 is similar to NGC 1042, where Ishields et al.1 j^Ooi) recently 
demonstrated the existence of a broad component in [N II], and that considering only the broad 
component moves the nucleus into the Seyfert /LINER regions in line-ratio diagnostic diagrams. 
Finally, the combination of a point-like hard source, soft extended circumnuclear emission, and an 
Fe Ka line with a large EW, is one often seen in type 2 AGNs. 

There is undeniably strong star formation occurring in the nucleus of NGC 4102, and we 
cannot rule out the extremely large SFR implied if all of the observed emission is imputed to star 
formation alone. Nevertheless, the evidence is strong that there is an AGN in NGC 4102. We 
conclude that NGC 4102 is another example where an AGN and strong star formation co-exist at 
the nucleus. 



3.4. NGC 4647 



NGC 4647 is a galaxy of ty pe Sc, its nucleus classified as H II in IHo et al.l (jl997al ). at a 
distance of 16.8 Mpc (|Tu11vIIi988| ). The Chandra observation is of the elhptical galaxy NGC 4649, 
and NGC 4647 is on the chip, at an off- axis angle of ~ 2.5'. Two factors complicate the detection 
of the nucleus of NGC 4647. First, the nucleus lies within the extended emission from NGC 4649. 
Second, the nucleus falls on a node boundary of the CCD. Wavdetect at its default filter for source 
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significance (~ 1 false source per 10^ pixels) does not detect the nucleus. However, once the elliptical 
galaxy is modeled and subtracted out of the image, there is a positive residual at the location of the 
nucleus of NGC 4647 (see Fig. [5]). To extract source counts we used a circular region centered on 
the centroid of the residual and with radius 2.3", which is approximately the 95% encircled-energy 
radius at 1.5 keV at that position. To extract background counts we located another region on 
the same node boundary that was at the same distance from the center of NGC 4649 as was the 
source circle. The source region has 11 counts after background subtraction, but this number 
necessarily has a large uncertainty. We take the radius of the source circle, 2.3", as the Chandra 
positional uncertainty of this source. 



The nucleus may also have been detected in x-rays by XMM-Newton (jRandall et al.ll2006l ). 
While the source was detected with S/N = 11, the positional uncertainty was 3". The XMM- 
Newton source is soft, similar to the Chandra residual. The XMM-Newton and Chandra source 
positions differ by 5.2". There is a nuclear 2MASS point source [Kg 12.3) at a distance of 1.1" 
from the Chandra position, but it is embedded in diffuse emission and the flux from the point 
source is poorly constrained. In the radio, there is a 5fT upper limit to the nuclear emission of 0.5 
mJy at 5 GHz (jUlvestad &: Hdl2002l ). Ther e is an older report of a radio de tection of the nucleus 
with a flux density of 16 mJy at 1.4 GHz (j Willis et al.lll976l : iKotanvil 11980 ). but with positional 
uncertainty {Aa,A6) = (0.16s, 11.2"). 



3.4-1- The nature of the nuclear emission 

The data are inconclusive at this time as to whether there is an AGN in NGC 4647. The 
faintness of the source and the positional uncertainties in the existing observations prevent a firm 
identification and characterization of the source. 



3.5. NGC 4713 



This galaxy is of type Sd, at a distance of 1 7.9 Mpc (lTullvlll988l l. Its nucleus was classified as 
T2 (transition object with type 2 spectrum) bv lHo et al.l (|l997al ). The nucleus is clearly detected 
by Chandra., but with only ten counts. It is a very so f t sour ce, with nine of the ten counts below 2.5 
keV. This object was also analyzed bv iDudik et al.l (|2005l ). but since they were looking for hard- 
band sources, this object was not counted as a detection. The observed count rate corresponds to a 
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flux of (5-10) xl0~^^ erg cm~^ s~^ in the 0.3-8 keV band, or luminosity (3-5) x 10^^ erg s~^ for our 
assumed distance. The 2-10 keV flux depends strongly on the assumed model, from /x ~ 1 x 10 
erg cm~'^ s~^ for a thermal bremsstrahlung model with kT = 0.3 keV, to fx 
s~^ for a power-law with T = 2 and Galactic absorption. 



5 X 10 erg cm ^ 



NGC 4713 was observed by HST in December 2006 (Martini et al. 2008, in preparation). 
The nucleus is resolved; thus only a nuclear star cluster or star forming region, and no AGN, is 
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detected. The cluster is ~ 0.40" in diameter in the image, corresponding to a physical size of ~ 35 
pc. Within this aperture the flux density is f\ = 5.34 x 10~^^ erg cm~^ s^^ which implies 

uL^ = 1.2 X lO""^ erg s"^ The nucleus is also in the 2MASS PSC, with J, H, Kg luminosities 
vLu ~ 10^^ erg s~^. However, the nucleus is embedded in diffuse extended emission and therefore 
the reported magnitudes may not be accurate estimates of the nuclear emissio n. There is an upper 
limit of 1.1 mJy to the nuclear radio emission at 15 GHz (jNagar et al.ll2005l ). and a similar limit, 
1.0 mJy, to the emission at 1.4 GHz from the FIRST survey. 



3.5.1. The nature of the nuclear emission 
As in the case of NGC 4647, the data are inconclusive regarding the presence of an AGN. The 



data ar e consistent with the AGN hypothesis: the "transition object" classification by iHo et al 



(|l997al ) implies there may be an AGN component in the optical spectrum; the nucleus is without 
doubt an x-ray source, though it is not possible to distinguish between emission from circumnuclear 
gas and the nucleus proper in the current Chandra imaging. The hardness ratio (with large uncer- 
tainty) is consistent with a F = 2 power-law. NGC 4713 may be similar to NGC 3184 in having a 
low-luminosity AGN inside a nuclear star cluster. 



3.6. NGC 5457 (M 101) 



NGC 5457 is a gala xy of type Scd at a distance of approxii nately 7 Mpc (IFreedman et al 



200ll : IStetson et al.l Il998l ) . The nucleus was classified as H II by IHo et al.l (jl997al ) . It has been 



observed multiple times by Chandra for a total observation time of about 1.1 Ms. Our analysis 
omits the shorter, and hence low signal-to-noise, exposures. The observations included here are 
listed in Table [3l The total usable exposure time is 695 ks. Chandra clearly resolves two sources in 
the nuclear region (Fig. [6]), whic h we label N and S . The northern source, N, is the nucleus, while 
the southern, S, is a star cluster ( Pence et ahlbooJ ). Table [3] shows the counts and hardness ratios 
of the two sources in each of the observations. Source N (the nucleus) varies in brightness by about 
a factor of 9 over the course of about 8 months (see Tableland Fig. [7]). There is no significant 
change in hardness ratio. In our analysis and discussion below we consider only the nucleus (source 
N). 

Three nuclear x-ray spectra were extracted: one ("merged") from the event list obtained by 
merging all observations done in 2004, and fit using the instrument response from ObsID 5339; 
one ( "high state" ) from merging two observations where the source had a high count rate (ObsIDs 
4736 and 6152) and fit using the response from ObsID 4736; and one ("low state") from merging 
five observations where the count rate was low (ObsIDs 5300, 5309, 4732, 5322, 5323), fit using the 
response from ObsID 4732. The spectra are shown in Figs. [8] and [9l fit models and parameters are 
shown in TableHl The low state spectrum can be fit with an absorbed power-law (x^/dof = 18.7/30) 
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but the fit is improved slightly with the addition of a plasma component (xsmekal) at temperature 
kT = 0.3 keV (Ax^ = —4.5 for two fewer d-o-f). The best-fit power-law slope T = 1.7ib0.5 is typical 
of unabsorbed AGN, and best-fit intrinsic absorption is consistent with zero. The unabsorbed 0.3-8 
keV luminosity is ~ 3 x 10^^ erg s~^. The high state spectrum can also be fit by an absorbed power- 
law. The best-fit slope is T = 2.2^q3, steeper but consistent with the low state slope within the 
uncertainties. The differences from the low state spectrum are, first, that the fit requires intrinsic 
absorption {N^ ~ 10^^ cm~^), and second, that there is no evidence for the MEKAL component. 
If a plasma component exists in the high state its flux falls below that of the power-law component. 
The unabsorbed 0.3-8 keV luminosity in the high state is ~ 3 x 10^^ erg s~^. Figure [TU] shows the 
high and low state spectra over-plotted. The spectra are identical below 1 keV, with the entire flux 
difference arising from the power-law normalization between 1 and 7 keV, and a possible "shoulder" 
between 1 and 2 keV in the high state. A fit where the power-law slope is forced to be identical in 
the low and high states is of similar statistical significance and produces best-fit values similar to 
the separate fits. This fit is shown in Table H] in the rows labeled "sim. low" and "sim. high" . The 
merged spectrum is similar to the high state one in that an absorbed power-law provides a good 
fit without the need for additional components. 

We analyzed archival HST WFPC2 images of NGC 5457, a 2400 s exposure using the F336W 
filter (referred to as U band below; nucleus on WF3) and a 1600 s exposure using the F547M filter 
(referred to as V band below; nucleus on PCI). The nucleus and the star cluster are detected in 
both images. In 2-pixel radius apertures, the U band nuclear fiux is 6.3 x 10^^'^ erg cm^^ s^^ 
A^^ and the V band flux is 8.3 x 10^^"^ erg cm~^ s~^ A~^, corresponding to m-pssew = 19-3 and 
= 19.0 in the STMAG system. The 2MASS PSC contains a source at the position of the 
nucleus (x-ray source N) but none at the location of source S. The given magnitudes J = 13.1, 
H = 12. 5, and K,^ = 11. 8 all correspond to uL^, « 10^^ erg s~^. The nucleus was not detected by 



FIRST. iHeckmanI (jl98d ) gives an upper limit to the nuclear luminosity at 6 cm corresponding to 



vLi, < 6 X 10^^ erg s ^ 



McElroyI ( 19951 ) reports cr* ~ 78 km s""*^ which implies log (M./Mq) « 2 x 10^. The corre- 
sponding Eddington luminosity is -LEdd ~ 3 x 10'^^ erg s^^. Assuming the source is an AGN, and 
assuming bolometric correction factors of ~ 10 and ~ 1 in the x-ray and IR, respectively, implies 
that the bolometric luminosity is in the range erg s ^, or that L/LEdd ~ 10-^-10-3. 



3.6.1. The nature of the nuclear emission 

The nuclear source shows several properties typical of an AGN. First, variability: The source 
varies by a factor of nine in eight months in 2004, and also varied between 2000 and 2004. Second, 
the x-ray spectrum: The spectrum is an absorbed power-law with slope ~ 2. In particular, the 
spectrum cannot be fit by a thermal component alone — a power-law is necessary. The low state 
spectrum is reminiscent of highly obscured AGN where an unabsorbed but diminished spectrum 
is seen via scattering. Third, x-ray colors: The ratio of 0.3-2 keV, 2-5 keV, and 5-8 keV counts 
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puts this object into the Compton-thick AGN part of the Levenson color-color diagram (Fig. 9 in 



Levenson et al.ll2006l ). In the context of the latter two points it is worth noting that if the source 
really is a highly absorbed AGN then energy emitted by the AGN may show up as re-processed 
radiation in the infrared. The existence of an infrared source with luminosity ~ 10^^ erg s~^ is 
consistent with this picture, though it does not argue for the presence of an AGN to exclusion of 
other types of sources. 

The x-ray properties of the nucleus are consistent with both the XRB and AGN hypotheses, but 
an AGN is not ruled out. The inferred x-ray luminosity (10'^''"^^ erg s~^) is in the range seen from 



XRBs, and the observed x-ray colors put this source in the LMXB region of the iPrestwich et al 



(j2003l ) color-color diagram. In favor of the XRB scenario is that the variability of the source is 
similar to the spectral state changes in XRBs. The low state may be interpreted as being the XRB 
state known as "hard" or "low/hard", and the best-fit power law slope, F ~ 1.7 is the value seen 
in XRBs in this state. The high state would then be the XRB state known as the "very high" 
or "steep power law" state. The best-fit power-law slope F ~ 2.2 in this case is less steep than is 
usually seen in XRBs in this state (F > 2.4) but the steeper value is included in the 90% confidence 
range. However, it must be kept in mind first, that at the quality of the spectra analyzed here the 
power law slopes are consistent with being iden tical in the two stat es, and second, that AGNs can 



show the same state behavior (e.g. IH 0419-577. iPounds et al.ll2004l ). In XRBs, the steep power-law 



state is associated with the presence of quasi-periodic oscillations in the x-ray emission, while the 
hard state is associated with the presence of a radio jet. In principle the presence of these features 
could provide additional evidence supporting the XRB hypothesis, but it is not currently feasible 
to detect them in XRBs at the distance of NGC 5457. 

Thus, there are again two possible scenarios, as in the case of NGC 3184. The x-ray source 
could be an HMXB in a super-star cluster. The star cluster would dominate the optical and IR 
emission. The large variation observed in the x-ray flux rules out the source's being more than one 
HMXB, as otherwise they would have to be varying in concert. The actual amount of obscuration is 
important, however, for the plausibility of the HMXB hypothesis since the inferred x-ray luminosity 
(~ 10^^ erg s~^ in the high state) is already at the high end of the range of XRB luminosities and 
there is not much room for a significantly higher absorption-corrected intrinsic luminosity. The 
alternative scenario is an AGN together with a nuclear star cluster. In AGNs both the intrinsic 
luminosity and the amount of obscuration are known to vary. The low state x-ray spectrum may 
be explained as a truly under-luminous (10^^ erg s~^) unabsorbed AGN, or as an AGN where the 
obscuration is so high that only scattered light, and no direct emission, is seen . The fact that the 



x-ray colors are similar to those of Compton-thick AGN ([Levenson et al.ll2006l ) supports the latter 
view, but we also note the absence of the 6.4 keV Fe Ka line in the x-ray spectrum that is often 
present in the reflected component. Though both the HMXB and AGN hypotheses are possible, 
on balance the AGN appears to be the more plausible one. 
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Discussion 



The motivation for this paper was to evaluate the feasibihty of detecting low-mass SMBHs 
in late-type spiral galaxies, which may still be accreting at the current epoch and if so should be 
detectable in x-rays. The six galaxies studied in this paper are not all late-type, but span the range 
Sa-Sd. NGC 3169 and NGC 4102 were regarded as low luminosity AGNs. None of the remaining 
four nuclei, however, was known to have an accreting SMBH, of any mass. NGC 3169 and NGC 
4102 are of type Sa and Sb, which have massive bulges and are expected to have massive SMBHs. 
For galaxies of types Scd and Sd, on the other hand, the lack of a luminous AGN could mean either 
that there is no SMBH or that the mass of the SMBH is low. As such, the observations studied 
here examine both aspects of a search for accreting low-mass SMBH: First, are these objects really 
detectable, given that the accretion rate is expected to be low? Second, are any sources detected 
in the very latest type spiral galaxies that have small or no bulges? 

We first note that of the six galaxies presented here, all six show nuclear x-ray sources. This 
implies that it is a very common occurrence. In a survey of late-type spiral galaxies such as 
our ongoing Chandra survey, therefore, the predominant concern is not going to be detection 
efficiency, but rather identification of the AGNs among the detected sources. Given that the sample 
presented in this paper consists of only six galaxies, we do not draw statistical conclusions here of the 
prevalence of very low-luminosity AGNs in nearby galaxies. But we note that, as shown in §§3.1-3.6, 
of the six nuclear x-ray sources, NGC 3169 is almost certainly an AGN, and NGC 4102, NGC 3184, 
and NGC 5457 , have very strong, though not conclusive, arguments in favor of their being AGNs. 
The two remaining galaxies, NGC 4713 and NGC 4647, are ambiguous but AGNs are not ruled out. 
We discuss below the issues such surveys will face when attempting to identify the nature of the 
detected sources. The diagnostic tools t raditionally used to distingui sh AGNs from non-AGNs (e.g. 



optical line ratios, iBaldwin et al. 



1981 



Veilleux &: Qsterbrocldll987l ) were developed in the course 



of studying luminous AGN. Dilution of the AGN emission by host galaxy light was not a serious 
problem and observations with low spatial resolution (several arcseconds) sufficed. In the study 
of AGN that are either intrinsically less luminous or are heavily obscured, howev er, host galaxy 



light becomes incre asingly problematic, and surveys relying on optical spectra (e.g. iHo et al 



1995 



Greene fc Hdl2004l ) require careful subtraction of the starlight using galactic spectral templates. In 
the weakest AGNs, however, signs of AGN emission may not be detected by the usual diagnostics. 
This problem will persist until optical observations with angular resolution of 1-10 mas become 
possible so that host galaxy light can be effectively excluded, though it can be mitigated by using 
regions of the spectrum where host galaxy emission is negligible, for example very high energy x-rays 
(tens to hundreds of keV). The detection of a compact radio source unresolved at milliarcsecond- 
scales, especially if the source is accompanied by jets, would also unambi guously ident i fy the nucleus 
as an AGN. This has been the motivation for radio surveys like that of lNagar et al.l (|2002l ). Some 
AGNs obsc ured in the optical and UV may be de t ectab le using infrared emission line strengths and 
ratios fe.g. bale etlHlbood : Isatvaual et ahlbooj boosl l. 



While the methods listed above allow the unambiguous identification of AGNs, observations 
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often do not have the angular resolution or sensitivity to distinguish AGN and non-AGN flux. 
Other sources of radiation in the vicinity of the AGN are, for example, plasma photoionized by 
the AGN itself, or a nuclear star cluster. The targeted AGNs have very low luminosity and thus 
even moderate amounts of obscuration may cause a significant decrement in the observed flux. In 
most cases, therefore, the AGN contribution should not be expected to dominate the total observed 
flux. Consequently, identifying these AGNs requires a different approach than what can be used 
in the case of the more luminous AGNs (Seyferts and QSOs). AGNs can be identified using x- 
ray observations (e.g. with Chandra and XMM-Newton) solely, but only if they are point sources 
whose inferred luminosities are greater than ~ 10"^^ erg s~^. Below that value, AGNs can become 
indistinguishable from ULXs and XRBs in x-rays. ULXs can have luminosities of a few times 



-|^q40 gj.g g 1 ^g_g_ Igoria et al. 2007 ) and have x-ray spectra that look similar to AGN spectra. It 



has been suggested (jStobbart et al.ll2006l ) that ULX spectra show a break at ~ 5 keV. AGNs are 
not known to show this break. In high quality spectra with a large number of counts it may be 
possible to exploit this difference to separate ULXs and AGNs. XRBs have power law spectra with 
r ~ 2, can show an Fe Ka emission line, and emit hard x-rays, all characteristics of AGNs as 
well. Additionally, even with Chandra^s angular resolution, the physical space probed ranges from 
the central ~ 10 to 100 pc of the galaxy. The existence of one or more XRBs within that region 
would be unsurprising. However, the fact that the inferred luminosities are as high as 10^^ erg s^^ 
severely constrains the expected number of XRBs. For example, fo r NGC 5457 , one of the four 



galaxies presented here that are not confirmed AGNs, iPence et al.l (|200ll ) provide a log N-log S 
relation as well as the surface density of x-ray point sources as a function of radius (their Figs. 3 
and 4). Approximately 12.5% of the sources have luminosities exceeding 10^^ erg s~^. The surface 
density of sources in the innermost 0.5' is ~ 4.75 arcmin"^. Therefore we may expect ~0.6 sources 
arcmin~^ above the luminosity cut-off in the central 0.5', or 4 x 10^'^ such sources within the 
Chandra source circle of radius 2.3" that has been used here. NGC 5457 is of type Scd, and can 
be taken to be representative of the other three galaxies, NGC 4647, NGC 3184, and NGC 4713 
which are types Sc, Scd, and Sd, respectively. Thus invoking XRBs and ULXs alone to explain 
the x-ray observations leads to physically implausible conditions, such as requiring that most, or 
all, quiescent, non-starburst spiral galaxies have a ULX or 100 XRBs in the central 0.5" — 2". 
Nevertheless, x-ray observations by themselves will usually be inadequate to distinguish AGNs from 
non-AGNs in any particular instance. 

Information from other wavebands is thus crucial for the identification process. For the six 
galaxies in this paper multi- wavelength photometry is summarized in Tables [5] and El But here too, 
traditional methods of identifying AGNs using fiux ratios such as aox, otxx, and fx/fn-, must be 
used with caution. The low luminosities of the AGNs mean that the emission in the two bands 
being compared may not be from the same object. For instance, for an obscured AGN surrounded 
by a nuclear star cluster, the observed x-ray fiux may be from the AGN but the optical fiux may 
be dominated by the cluster. The existence of an AGN must instead be inferred by consistency 
and plausibility checks considering as much of the spectral energy distribution as possible, and the 
goal is the rejection of the hypothesis that all of the observed properties can be explained without 
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requiring the presence of an AGN. NGC 3184 provides a good example where different modes of 
observing, imaging and spectroscopy, in two wavebands, x-ray and IR, together make a compelhng 
argument for the presence of an AGN where each observation individuahy is ambiguous. 



Nuclear star clust ers are fairly common in spiral galaxies (e.g. iBoker et al.ll2002l : IWalcher et al 



20051 : ISeth et al.ll2008l ). A cluster poses two main problems. First, it makes the presence of XRBs 
more likely. Second, if the cluster is young and contains many O and B stars, it may o verwhelm 
AGN emission in the UV in addition to the optical (e.g. NGC 4303, IColina et al.l 120021'). It may 



be po ssible in some cases, as for NGC 1042 (jShields et al.ll2008l ) and NGC 4102 (| Gonial ves et al 



19991 ). to attempt a separation of the cluster and AGN components in the optical spectrum. In 
addition, stellar spectra, even for late-type stars, fall faster towards the infrared than the power- 
laws typical of AGNs. The mid-infrared colors of AGNs, therefore, tend to be redder t han those of 



stella r populations and this color difference can be used to infer the presence of an AGN (jStern et al 
200rI ). 



In addition to the spectral energy distribution, source variability can be a discriminant, as 
AGNs are known to vary in all wavelengths, whereas a star cluster, say, would not. Conversely, if 
the variation is periodic it would rule out an AGN and argue for an XRB instead. 

A survey of the type discussed here finds AGNs and strong AGN candidates, but does not mea- 
sure the masses of the SMBHs in those AGNs. Measurement of the mass of one of these SMBHs 
will be a difficult endeavor, since the sphere of influence of the black hole cannot be resolved with 
current technology and resources. None of the six objects studied in this paper shows broad op- 
tical emission lines whose widths could be used to estimate the BH mass, and this is likely to be 
typical behavior. Spectropolarimetry may uncover broad lines in polarized light in some of the 
AGNs. Estimates of the SMBH mass may be made by using known scaling relationships, with the 
caveat that the correlations are all based on SMBHs two or more orders of magnitude more massive 
than the ones expected to be found by the survey. The least indirect method is ai i application of 



the o bserved correlation between the x-ray power law slope and Eddington ratio (IWilliams et al 



2004I ). The Eddington ratio and the luminosity in turn provide an es timate of the mass of the 



SMBH. Other relatio nships that can be used are: (a) the M,-a relation (jFerrarese &: Merritt 



200C 



Gebhardt et al.ll200d ). but this method becomes more and more uncertai n as the bulge itself be- 
comes ill-defined in the latest-type spirals; (b) the M, -Lhu]^c relationship (jKormendv &: Richstone 



1995 



McLure fc Dunlop 



2001 



Marconi &: HuntI |2003| ) , which h as more scatter and also faces th e 
problem of the definition of the bulge; (c) the iVf,— fcirc 

relation (|Ferraresell2002l : iBaes et al.ll2003l ). 



which ha s the advantage that it does not require the presence of a well-defined bulge; (d) the M,-C 
relation (iGraham et al.l l200ll ) , where C is the concentrat ion of light. There is also a reported re - 
lationship between black hole mass and core radio power (ISnellen et al.ll2003l : iMcLure et al.ll2004l ). 
but this relationship is not as well established as the others, and is based on observations of ellip- 
tical galaxies, so its applicability to the spiral galaxies here is uncertain. Overall, there is unlikely 
to be one standard method of measurement that can be applied to these galaxies, but one or more 
of the above methods may provide useful estimates of or limits to the masses of the SMBHs in 
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these AGNs. Mass measurements independent of the scaHng relationships are possible if an object 
turns out to have broad emission lines, like NGC 4395, in which case line widths or reverberation 
mapping may be used, or if it has maser emission, like NGC 4258, in which case gas dynamics can 
be used. Mass measurement in other cases will have to await mas-scale angular resolution in bands 
other than radio to resolve the sphere of influence of these black holes. 

Of the six galaxies here, NGC 3169, NGC 4102, and NGC 5457 have measurements of either 
the stellar velocity dispersion or the luminosity of the bulge, thus allowing an estimate of their 
SMBH masses (assuming here that the source in M 101 is an AGN). The scatter in the M,-a 
and Af,-Lbuige relations, together with the uncertainty in the observed flux and the bolometric 
correction, result in uncertainties of about an order of magnitude in the inferred Eddington ratio, 
but all three objects have L/Lkhh ^ 10~^- T his is in the range seen in low-luminosity AGNs (e.g. 
L ~ 10~^LEdd for M 81; lYoung et al.l 120071 ). and much higher than the L ~ 10~^LEdd seen in 



truly quiescent SMBHs. These observations indicate that there is indeed a population of accreting 
SMBHs in nearby spiral galaxies that do not show optical signs of activity but can be uncovered by 
looking for their x-ray emission. Such a population will answer the question of whether the bulge 
is the dominant component that determines the existence, and mass, of a nucle ar SMBH. The 



disco veries of AGNs in the Sd galaxies NGC 4395 (|Ho et al.lll997bl ') and NGC 3621 (jSatvaDal et al 



20071 ). and the strong evidence for AGNs in the Scd galaxies NGC 3184 and NGC 5457 suggest 
it is not, at least as far as existence is concerned. Among the SMBHs discovered in the latest- 
type spirals (with small or no bulges) and in the lowest-mass galaxies should be a population of 
SMBHs with masses less than 10^ Mq, enabling a systematic study of the low-mass end of the local 
supermassive black hole mass function. 
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1 2 3 4 5 6 7 8 




Energy (keV) 

Fig. 1. — Spectrum of NGC 3169. Squares show the data binned by 20 PHA channels. The soUd 
hue shows the best-fit model, an absorbed power-law with Nh ~ 10^^ cm~^, and T = 2.6. The 
lower panel shows the residuals from the fit. 



Table 1. Targets and observation parameters 



Target 


Morph. 


Nucleus 


Coordinates (J2000) 


Dist.*^ 


Scale 


Obs. Date 


ObsID 


Exp. 




Type 


Type'^ 


RA 


Dec 


(Mpc) 


(pc/") 






(ks) 


NGC 3169 


Sa 


L2 


10 14 15.0 


+03 27 58 


19.7 


96 


2001 May 2 


1614 


2.0 


NGC 3184 


Sod 


hi: 


10 18 17.0 


+41 25 28 


8.7 


42 


2000 Jan 8 


804 


39.8 
















2000 Feb 3 


1520 


21.3 


NGC 4102 


Sb 


hi: 


12 06 23.1 


+52 42 39 


17.0 


82 


2003 Apr 30 


4014 


4.9 


NGC 4647 


So 


HII 


12 43 32.3 


+11 34 55 


16.8 


81 


2000 Apr 20 


785 


36.9 


NGC 4713 


Sd 


T2 


12 49 57.9 


+05 18 41 


17.9 


87 


2003 Jan 28 


4019 


4.9 


NGC 5457^= 


Scd 


HII 


14 03 12.6 


+54 20 57 


7.2 


35 


C 


C 


750 



'^From lHo et al.l |l997al ). L2: Type 2 LINER, T2: Type 2 Transition object. 
'^FromlTuUd (|l988h except from Istetson et al.l (Il998h for NGC 5457. 

'^NGC 5457 (M 101) was observed multiple times. The total observation time analyzed here is given in this 
table. The individual observations are listed in Table [3] 
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1 2 3 

Energy (keV) 



Fig. 2. — Left: Nucleus of NGC 3184, which appears to consist of two components, labeled N and S 
in the image. North is up and east to the left in the image. About 117 and 36 counts were detected 
in the southern and northern components, respectively. The image is 15" on a side. The black bar 
in the lower left represents a projected distance of 4", corresponding to ~ 170 pc. Right: Spectrum 
of the southern component of the nucleus of NGC 3184, together with an absorbed power-law 
model. Data points were binned so that there were at least 5 counts in each bin. The model shown 
above (red line) has A'^^ = 3.0 x 10^^ cm^'^ and F = 2.6. The bottom panel shows the residuals. 
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Fig. 3. — Nucleus of NGC 4102 with the regions used to extract counts overlaid on a hardness ratio 
map. North is up and east to the left in the image. The lightest colored pixels have HR = — 1; 
the darkest pixels have HR = +1. The circle, 2.9 pixels in radius, was used for the core. The 
partial annulus is the region used to extract counts from the extended emission. The compact core 
has harder emission than the extended component. The bar on the lower left represents 2", or a 
projected distance of ~ 160 pc. The image is 10" on a side. 
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Table 2. X-ray Measurements 



Target Counts Bkg/Src HR' 

Broad Hard Soft Area Ratio 





Src 


Bkg 


Src 


Bkg 


Src 


Bkg 








NGC 3169 


159 


23 


148 


5 


11 


18 


21.4 


+0.86 


+0.05 
-0.03 


NGC 3184 N 


36 


95 


4 


33 


32 


62 


92 


-0.75 


+0.08 
-0.13 


S 


117 


95 


13 


33 


104 


62 


60 


-0.77 


+0.05 
-0.07 


NGC 4102 all 


354 


48 


78 


8 


276 


40 


6.8 


-0.55 


+0.04 
-0.05 


core 


171 


48 


68 


8 


103 


40 


52.5 


-0.20 


+0.08 
-0.07 


ext 


115 


48 


6 


8 


109 


40 


26.1 


-0.88 


+0.03 
-0.06 


NGC 4647 


15 


38 


1 


15 


14 


23 


10.1 


-0.80 


+0.04 
-0.20 


NGC 4713 


10 


4 


1 





9 


4 


21.4 


-0.69 


+0.09 
-0.25 


NGC 5457 


314 


256 


23 


45 


291 


211 


20.8 


-0.86 ± 0.03 



'^Hardness ratio, HR = (H— S)/(H+S), where H and S are the 
counts in the hard and soft bands respectively, calculated using the 
method described in iPark et al.l (|2006l ) . The tool used is available at 



http : //hea-www . harvard . edu/AstroStat/BEHR/[ 



Table 3. NGC 5457 X-ray Measurements 



ObsID 


Obs Date 


Exp Time 




North source 




South source 








Counts 


Rate 


HR 


Counts 


Rate 


HR 






(ks) 




[KS ) 






[KS ) 




934 


2000 Mar 26 


97.8 


262 


2.68 


n 0-1+0.03 
-U.»i_o.o4 


264 


2.70 


-0.70 ± 0.04 


4731^ 


2004 Jan 19 


56.2 


92 


1.64 


-0.64in^ 


136 


2.41 


-0.58t°:°? 


5300 


2004 Mar 7 


52.1 


36 


0.69 


-0.75lH^ 


117 


2.24 




5309 


2004 Mar 14 


70.8 


38 


0.54 


-0 62+°-^° 


158 


2.24 


-0.641°:°^ 


4732 


2004 Mar 19 


69.8 


56 


0.80 


-0 77+0-0''' 
'-0.09 


151 


2.17 


-0.581°:°? 


5322 


2004 May 3 


64.7 


54 


0.83 


_0 70+0.06 
'-'•'^-0.10 


144 


2.22 


-0.69 ± 0.06 


5323 


2004 May 9 


42.4 


42 


0.99 


-0 64+°-°^ 

^•"^-0.14 


120 


2.83 


^•-^^-o.os 


5338^ 


2004 Jul 6 


28.6 


61 


2.14 


-0 f^fi+0-09 


41 


1.45 




5339^ 


2004 Jul 7 


14.0 


18 


1.26 


-0 60+°-i3 

U.OU_Q 22 


27 


1.90 


_n 70+0.06 
^•'^-0.15 


5340 


2004 Jul 9 


54.4 


81 


1.48 


r, 79+0.06 

"-"•'^-o.os 


110 


2.02 


_n 74+0.04 


4734 


2004 Jul 11 


35.5 


31 


0.87 




63 


1.77 


-0 61+0-09 


6114 


2004 Sep 5 


38.2 


91 


2.38 


-0 60+°-°^ 

^•"'-'-0.09 


64 


1.68 


-0.641H? 


6115 


2004 Sep 8 


35.4 


96 


2.72 


n 78+0.05 


63 


1.79 


-0 70+0-07 

"-'-0.10 


4735 


2004 Sep 12 


28.8 


70 


2.43 


^•'"-0.10 


49 


1.71 


n 77+0.06 
' ' -0.11 


4736'^ 


2004 Nov 1 


77.4 


290 


3.75 


n 7Q+0.03 


122 


1.58 


-0.671°:°? 


6152*^ 


2004 Nov 7 


26.7 


133 


4.99 


n 7O+0.05 


51 


1.92 


-0.631°:°^ 



'^Source extended, or possibly image smeared. 
''Possible residual contamination from background flare. 



Table 4. X-Ray Spectral Fits 



A. dl t 


i-l-Si 1 oi'io 
vjr ojI ojL/ u 1*^ 






Spectral Fit Parameters 










Nh 




kT 


Nh 


r 


Refl.'^ 


Line 


EW 






(10^0 cm-2) 




(keV) 


/ 1 n21 — 2\ 






(keV) 


(keV) 




NGC 3169 


2.86 


ab(pl) 
ab(pl) 




00+46 
^^-36 

1161^° 


2.011:? 

2 6+^1 








1.1(26)^^ 
0.3(22) 


NGC 4102 core 


1.79 


ab(pl) 
abCrnj 




0+1.4 
0+0.4 


2.0 ±0.5 
1.8 ±0.4 


108lif 






7.4(27)'^ 
2.3(26)^^ 






abCrn+gj 




0+1.9 


rt, i^-UO ft 

^•"'-o.s 


129^1 


6.4 


3.3 


0.6(25) 






ab(rn+g) 




0+0.6 


9 9+0.6 


12911' 


6.4 


2.5 


1.6(25)'= 


ext 




ab(br) 


1 9+8.6 
-^•^-0.8 


1 9+2.8 
■■-•^-1.2 










0.5(18) 


NGC 5457 sep. low 


1.15 


ab(me+pl) 


0.31°:? 


0+1.8 


1.7 ±0.5 








0.5(28) 


sep. high 




ab(pl) 




■■-•^-1.5 


9 9+0.4 
^•^-0.3 








0.5(40) 


sim. low 




ab(ine+pl) 


0.31°:? 


42+1-9 
'-'■^^-0.42 


2 0+°-"^ 








0.5(28) 


sim. high 




ab(pl) 




1.1 ±0.6 


Z.U_Q 2 








0.5(41) 


merged 




ab(pl) 




0.6 it 0.4 


1.9 ±0.2 








0.6(70) 



^Model labels: ab=xswabs, photoelectric absorption; ga=xswabs with value frozen at Galactic column density towards 
this target; br=xsbremss, thermal bremsstrahlung; g=gaussld, one-dimensional Gaussian; me=xsmekal, thermal plasma; 
pl=powlawld, one-dimensional power law; rn=xspexrav, power-law reflected by neutral material. 

'^Reflection scaling factor. 

'^Cash statistic (cstat in Sherpa). 

'lUnconstrained by fit. 
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Table 5. Inferred nuclear luminosities 





Filter 


NGC 3169 


NGC 4102 


NGC 4647 


NGC 3184 


NGC 5457 


NGC 4713 




or Band 


(Sa) 


(Sb) 


(So) 


(Scd) 


(Scd) 


(Sd) 


X-ray 


0.3-8 keV 


41.7 


40.2 


39.0^ 


37.3^ 


37.5-38.5^ 


38.6^ 


UV & Opt. 


F300W 

F336W 
F547M 
F606W 








< 39.4 
39.5 


39.0 
39.5 


40.1 


IR 


Ks 


42.6 


42.8 


38.8 


40.9 


41.0 


41.0 


Radio 


15 GHz 
5 GHz 


37.2 








< 35.8 


< 36.8 




1.1 GHz 




38.0 


36. 9^' 


< 35.1 


< 31.9 


< 35.7 



Note. — Values are log(L/erg s^"*^) in the specified bandpass in the x-ray and log(z^L,y/erg s~^) where 
a single filter, wavelength, or frequency is given. Note that the values are derived from observations with 
varied instruments, PSFs, apertures, and signal-to-noise ratios, and the uncertainties can be as large as 
a factor of two. 

'"0.3-12 kcV. 

'^Large positional uncertainty makes it unclear whether the detected source is really the nucleus. 

'^Based on a small number of counts. A power-law spectrum was assumed. Intrinsic absorption is 
unknown. 

'^Variable source. 
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Fig. 4. — Spectra of NGC 4102. Data have been grouped to 5 counts per bin. Left: Spectrum of the 
extended component. The model is absorbed bremsstrahlung emission, with Nh = (l-2?;?.2) X lO^i 
cm~^ and kT = 1.2^'^Q'g keV. Right: Spectrum of the core. The model is a reflected power-law and 
a Gaussian line, with no absorption in excess of Galactic. The best-fit parameters are T ~ 2.3 it 0.6 
and reflection factor R « 130. The hue was fixed at 6.4 keV with FWHM = 0.3 keV. 



Table 6. NGC 3184 luminosity in Spitzer bands. 



Band 
(urn) 


Aperture 
Radius 


log(z^L^) ^ 
(erg s^-*^) 


3.6 


3" 


40.6 


4.5 


3" 


40.4 


5.8 


3" 


40.7 


8.0 


3" 


40.8 


24 


6" 


41.1 


70 


10" 


41.5 


160 


10" 


41.7 



'^Uncertainties are 10% in the 
3.6//m, 4.5;um, 5.8/im, 8.0/im, 
and 160/xm bands; 4% at 24/im; 
12% at 70/im. 
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Fig. 5. — On the left are shown the residuals in the Chandra image after the elliptical galaxy NGC 
4649 is modeled and subtracted. The residual image has been smoothed by convolving it with a 
Gaussian with a = 5 pixels. The black circle, whose radius is 2.3" or roughly the 95% encircled- 
energy radius at 1.5 keV at that position, shows the probable nucleus of NGC 4647. On the right, 
a circle (white) of the same radius, whose center has the same celestial coordinates as the one on 
the left, is superposed on a DSS image of NGC 4647. Both images are shown on the same scale 
and have north up and east to the left. The bar in the lower left corresponds to 10", or a projected 
distance of ~ 800 pc. 
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Fig. 6." Nucleus of NGC 5457 (M 101). North is up and East to the left in the image. The 
nucleus is resolved into two sources, marked N and S in the image. The northern source N is the 
candidate active nucleus. The southern source S is a known star cluster (jPence et al.ll200ll ). The 
black bar in the lower left represents a projected distance of 2", or ~ 70 pc. The image is 8"on a 
side. 
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Fig. 7. — Variability and hardness ratio of the two NGC 5457 sources. Only observations performed 
in 2004 are shown here. The nucleus (source N) varies by a factor of ~ 9 between March and 
November. Uncertainty in count rate was derived assuming y/n uncertainty in the counts. 
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Fig. 8. — Spectrum of NGC 5457 obtained from merging all 2004 observations. Data have been 
binned to 15 counts per bin. The line shows the best fit model, an absorbed power-law (F = 1.9) 
and thermal plasma (kT = 0.44 keV). The line-like feature at 4 keV may be due to statistical 
fluctuation. The lower panel shows the residuals. 




Energy (keV) Energy (keV) 

Fig. 9. — Spectrum of NGC 5457 in low and high states. Data have been binned to five counts per 
bin. Left: Low state spectrum. The best-fit model is a power-law plus MEKAL plasma with no 
absorption in excess of Galactic. Right: High state spectrum. The best-fit model is a power-law 
with both intrinsic and Galactic absorption but no plasma. 
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Fig. 10. — Spectrum of NGC 5457 in the high state (red) and the low state (green) over-plotted. 
The spectra are similar below 1 keV. Most of the difference is in the hard component. 
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